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Abstract 



p • We study complex-temperature properties of the uniform and staggered susceptibilities x 

and x^""^ of Ising model on the honeycomb lattice. From an analysis of low-temperature 
series expansions, we find evidence that x ^ind x^^^ both have divergent singularities at 
I the point z = —1 = Z£ (where z = e~'^^), with exponents 7^ = 7^,j = 5/2. The critical 
amplitudes at this singularity are calculated. Using exact results, we extract the behaviour 
of the magnetisation M and specific heat C at complex-temperature singularities. We find 
that, in addition to its zero at the physical critical point, M diverges at z = — 1 with 
exponent (3i = —1/4, vanishes continuously at 2; = ±i with exponent f3s = 3/8, and vanishes 
discontinuously elsewhere along the boundary of the complex-temperature ferromagnetic 
phase. C diverges at z = —1 with exponent a'^ = 2 and at t> = ±i/^/3 (where v = tanhi^) 
with exponent cte = 1, and diverges logarithmically at z = ±i. We find that the exponent 
relation «' -|- 2/3 + 7' = 2 is violated at 2; = — 1; the right-hand side is 4 rather than 2. The 
connections of these results with complex-temperature properties of the Ising model on the 
triangular lattice are discussed. 
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1 Introduction 



There are several reasons for studying the properties of statistical mechanical models with 
the temperature variable generalised to take on complex values. First, one can understand 
more deeply the behaviour of various thermodynamic quantities by seeing how they behave 
as analytic functions of complex temperature. Second, one can see how the physical phases 
of a given model generalise to regions in appropriate complex-temperature variables. Third, 
a knowledge of the complex-temperature singularities of quantities which have not been 
calculated exactly helps in the search for exact, closed- form expressions for these quantities. 
This applies, in particular, to the susceptibility of the 2D Ising model, which has, to this 
day, never been calculated, in contrast to the (zero-field) free energy, first calculated (for the 
square lattice) by Onsager and the spontaneous magnetisation, the expression for which 
(for the square lattice) was proposed by Onsager and first calculated by Yang p| . A fourth 
reason for the interest in complex-temperature singularities is that they can significantly 
influence the behaviour of a given quantity for physical values of the temperature. Indeed, 
early studies 0, ^, |^ of such complex-temperature singularities were motivated in part by 
the fact that when they occurred closer to the origin, in a certain low-temperature expansion 
variable, than the physical critical point, they precluded the application of the ratio test, 
then in common use, to determine the location of this critical point. The first work on 
natural boundaries of the free energy of the Ising model (on the square lattice) for complex 
temperature was in Ref. 0. 

Several years ago, we reported some results on complex-temperature singularities of the 
susceptibility and correlation length for the Ising model , including a discussion of complex- 
temperature symmetries and a proof (for the square lattice) that the (zero- field) susceptibility 
can have at most finite non-analyticities on the border of the complex-temperature extension 
of the symmetric phase, apart from its divergence at the physical critical point. Two recent 
papers have reported results on complex-temperature singularities for the Ising model on 
the square lattice p|, ||. In Ref. [|l^, the present authors extended the study by Guttmann 
1^ of complex-temperature properties of the Ising model on the triangular lattice. In the 
present paper we shall carry out a study of the complex-temperature properties of the Ising 
model on the honeycomb (= hexagonal) lattice. For completeness we note that an analysis 
of complex-temperature singularities in the 3D Ising model, after those of Refs. 0, Q , was 
given in Ref. |TT[| . 
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2 Complex- Temperature Extensions of Physical Phases 



Our notation follows that in our previous papers p, [T^, so we review it here only briefly. 
We consider the Ising model on the honeycomb lattice (coordination number g = 3) at 
a temperature T and external magnetic field H defined by the partition function Z = 
J2{(7,} with the Hamiltonian 

n = -J ^^^J-HY.ai (2.1) 

<ij> i 

where ai = ±1 are the Z2 spin variables on each site i of the lattice f3 = {kBT)~^, and J 
is the exchange constant. We use the standard notation K = (3 J, h = f3H, v = tanhK, 
z = e-^^ = (1 - v)/{l + v), and u = z . Another relevant variable is the elliptic modulus 
which one encounters in the broken-symmetry phases, 

_ 4.3/2(1 + , 2)1/2 

^< " {l-z)%l + z) ^^-^^ 
and its inverse, which occurs in expressions in the Z2-symmetric phase, 

^> = = (1_ ^2)3/2(1 +3^2)1/2 (2-3) 

We note the symmetries 

K ^ -K =^ {v ^ -V , z^l/z, u^l/u, k^^-k^} (2.4) 

where k^ = k^ or /c>. The reduced free energy per site is / = — /5F = lim^v^^oo A^^^"^ In Z 
in the thermodynamic limit, where Ng is the number of sites on the lattice. The zero-field 
susceptibility is x = ^^Qjf ^ \h=o, where M{H) denotes the magnetisation. The staggered 
susceptibility is denoted x*^"^^- It is convenient to deal with the reduced quantities x = P^^X 
and x*-"^ = P'^X^""^- We recall that on a loose-packed lattice such as the honeycomb lattice, 
in the symmetric, paramagnetic (PM) phase, the uniform and staggered susceptibilities are 
simply related according to 

X^''\v) = xi-v) (2.5) 
Following Onsager's solution for f{K, h = 0) on the square lattice, the free energy was cal- 



culated for the honeycomb lattice in the papers of Ref. ||T2|. The spontaneous magnetisation 
for the honeycomb lattice was first given by Naya [Q. The critical coupling separating the 
symmetric, paramagnetic (PM) high-temperature phase from the phase with spontaneously 
broken Z2 symmetry and ferromagnetic (FM) long-range order is Kc = (1/4) In 3, so that 
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Vc = l/-\/3 = 0.577350... and Zc = 2 — ^/S = 0.267949... As usual for loose-packed lattices, 
the critical point separating the PM phase from the phase with antiferromagnetic (AFM) 
long-range order is = —Kc, or equivalently, v = —Vc, z = l/zc = 2 + ^/S. 

We begin by discussing the phase boundaries of the model as a function of complex 
temperature, i.e. the locus of points across which the free energy is non-analytic. As noted 
in Ref. 0, there is an infinite periodicity in complex K under certain shifts along the 
imaginary K axis as a consequence of the fact that the spin-spin interaction a^aj in 7i is an 
integer. In particular, there is an infinite repetition of phases as functions of complex K; 
these repeated phases are reduced to a single set by using the variables v or z owing to the 
symmetry relation K —>■ K + nm ^ {f — > f , z z}. The requisite complex extensions of 
the physical phases can be seen by using the exact expression for the free energy [|T^ , 

/ = ln2 + - 1^ j ^ \n[-[cosh\2K) + 1 - sm]i\2K)P{e,, 6,)] } (2.6) 



where 



P(^i,^2) = cos^i + cos^2 + cos(^i + ^2) (2.7) 

The boundaries of the complex-temperature phases are comprised of the locus of points where 
the argument of the logarithm in / vanishes.^] Expressed in terms of the low-temperature 
variable z, / is given by 

/ = |ir+l££^ln[(l + .)^{(l-2. + 6z^-2z^ + ;.^)-2z(l-^)V(^i,^^2)}] (2.8) 

(where g = 3). Evidently, the argument of the logarithm vanishes along the curve defined 
by the solutions to the equation 

{l-2z + Qz^ -2z^ + z^)-2z{l~ zfx = Q , < a; < 3 (2.9) 

where x = P{6i, 62). This curve is shown in Fig. 1(a). (Note that the curve contains the point 
z = —1 where the initial factor (1 + z^ vanishes.) Since eq. (|2.9|) has real coefficients, the 
solutions are either real or consist of complex conjugate pairs, which explains the reflection 
symmetry of the curve in Fig. 1(a) about the real axis in the z plane. Furthermore, under 
the transformation z 1/z, the left-hand side of eq. (p.9|) retains its form, up to an overall 
factor of z~^] consequently the locus of solutions given by the curve in Fig. 1(a) is invariant 
under this mapping, z 1/z. For x = 3, eq. (^I9|) has double roots at z = Zc, l/zc. As x 
decreases from 3 to —1, the pairs of complex solutions move along the curve, rejoining again 



^ The free energy is trivially infinite aX K = 00; since this is an isolated point and hence does not form 
part of a boundary separating phases, it will not be important here. 
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in two pairs of double roots at z = ±i for x = —1. Finally, as x decreases from —1 to —3/2, 
the solutions move outward from z = ±i along the unit circle; the leftward-moving roots 
join at z = —I while the rightward-moving roots terminate at the endpoints z = e^™!"^ . We 
shall denote the leftmost (£) of the real roots as z^ = —\. The points z = are singular 
points of the curve in the technical terminology of algebraic geometry specifically, they 
are multiple points of index 2, where two arcs of the curve cross each other (with an angle 
of vr/2). We denote these as Zs±-, respectively. The endpoints at z = e^*'^/^ are, of course, 
also singular points of the curve in the mathematical sense. The physical phases are: (i) 
FM, for < z < Zc, (ii) PM, hi z^ < z < 1; and (iii) AFM, for 1 < z < oo. The 
complex-temperature extensions of these are the regions marked in Fig. 1(a). Note that the 
complex-temperature phase boundaries in Fig. 1(a) are formally the same as those for the 
Ising model on the triangular lattice, in the v plane (shown in Fig. 1(c) in Ref. |0). This 
is a consequence of the geometric duality of the honeycomb and triangular lattices. Recall, 
however, that the actual phase structures, both for physical and complex temperature, are 
different; in particular, the model on the triangular lattice has no AFM phase, and the outer 
phase in the v plane (denoted O in Ref. |10|) has no overlap with any physical phase. 

The transformation from z to the elliptic modulus /c< given in eq. ( p.2|) maps the z 
plane to the image shown in Fig. 1(b). As noted in ( p.4|) , z and 1/z are mapped to A;< 
and — fc<, respectively. The point Zc is mapped to A;< = 1. The three points l/z^ and ±i 
are all mapped to A;< = —1. The point 2; = 0, as well as the points at complex infinity 
(i.e., z = limp^oope*^ for arbitrary real 6) are mapped to fc< = 0. In addition, given the 
factorisation (1 — z+z'^) = (z — e^^^^){z — e~^'^^^), it follows that the finite endpoints z = e^'^'^l'^ 
are mapped to fc< = 0. This latter property has the consequence that this elliptic modulus 
variable is not a useful one in which to re-express the low-temperature series expansion for 
X, in contrast to the case of the square lattice, where the analogous variable /c< played a 
very valuable role. The transformation given by (|2.2|) maps the complex-temperature FM 
and AFM phases formally onto the same phase, occupying the interior of the unit circle in 
the /c< plane; similarly, it maps the complex-temperature PM phase onto the exterior of this 
circle. The complex conjugate arcs of the unit circle from arg(z) = ±7r/3 to arg(z) = ±7r/2 
are mapped onto the same line segment extending along the negative real axis in the (fc<) 
plane from to —1, while the complex conjugate arcs of the unit circle from arg(2) = ±7r/2 
to z = —\ are mapped onto the negative real axis from A;< = — 1 to A;< = —00. As usual 
with a singular point, the image of the point z = —\ depends on the direction in which one 
approaches it in the z plane; for example, if one approaches this point along the negative 
real z axis (with a infinitesimal positive imaginary part so that with the branch cut for the 
z^l"^ factor being placed along the negative real axis, (—1)^/^ evaluates to — i), the image 
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/ = In 2 - ^ ln(l -v') + \ j\ j\ ln[l + - 2v'il - v')P{e,, 6,)] (2.10) 



approaches fc< = — ioo, and so forth for other directions. 

One may work out the phase boundaries for the honeycomb lattice in the v plane either 
by transforming the boundaries in the z plane using the bilinear conformal mapping v = 
(1 — z)/{l + z), OT directly by analysing the free energy expressed in terms of v, viz., 

(The second term is just the usual (g/2) ln(coshi^') term.) Aside from the trivial infinity 
at f = 1 {K = oo), the locus of singularities of / is given by the solutions to the equation 
l+3f 2f ^(1— t'^)x = for —3/2 < x < 3. This is shown in Fig. 1(c). This equation has real 
coefficients and is invariant under the transformation v —>■ —v; the respective consequences 
of these two properties are that the solution curve in Fig. 1(c) is invariant under reflection 
about the horizontal and vertical axes in the v plane. The vertical line segments extend 
from ±i/^/3 to ±ioo, respectively, and are the images, under the conformal mapping, of the 
portion of boundary in the z plane lying on the unit circle. The curve in Fig. 1(c) is formally 
the same as the complex-temperature phase boundary for the Ising model on the triangular 
lattice in the variable z, because of the duality noted above. 

Using the general fact that the high-temperature and (for discrete spin models such as 
the Ising model) the low-temperature expansions have finite radii of convergence, we can use 
standard analytic continuation arguments to establish that in addition to the free energy and 
its derivatives, also the magnetisation and susceptibility are analytic functions within each 
of the complex-temperature phases. This defines these functions as analytic functions of the 
respective complex variable {K, v, or z, etc.). Of course, these functions are, in general, 
complex away from the physical line — oo < < oo. 

Our definition of singular forms of a function at a complex-temperature singular point was 
given in Ref. 0, |10|]. Note, in particular, that whereas a physical critical point can only be 
approached from two different phases, high- or low-temperature, some complex-temperature 
singular points may be approached from more than two phases. 

The low-temperature series for the staggered susceptibility is expressed in terms of the 
variable y = 1/z, and for our analysis of this series, we observe that the complex-temperature 
phase diagram in the y plane has the same phase boundaries as those in Fig. 1(a), owing to 
the invariance of this boundary under z ^ 1/z. The phases are, of course, inverted, so that 
the innermost phase is AFM, to its right, PM, and, in the outer region, FM. 

Finally, because of the star-triangle relations connecting the Ising model on the triangular 



and honeycomb lattices, the following exact relations hold [pO| : 



Xt[u) 



Ir . (a). 



xJ^') + x'HJi^') (2.11) 
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where 

■ 



z' 



u = (2.12) 

1 -.1 1^/2 ^ I 



and 



1 



XM = ^ Xnci^) + xSMJ = ^[Xh.iv) + x^.(-^^)J (2.13) 



1 



2 

where 

= ^^V^ (2-14) 
1 — w + 

3 Complex- Temperature Behaviour of the Specific Heat 
3.1 General 



From the exact expression (|2.8| ) for the free energy / we calculate the specific heat in the 
low-temperature phase as 

i.-^T^-^r- ' ~ ^2(2: + z') + 28(z^ + z') - 20{z' + z') + 18z'] ^ 

-^^^Eik^) (3.1.1) 
7r(l — 2; + z^) 

where K{k) = ^'"^ dQ\\ - P sin^ ^]-V2 and E{k) = /o"^^ rf^[l - sin^ ^] V2 are the complete 
elliptic integrals of the first and second kinds, respectively, which depend on the (square of 
the) elliptic modulus k, and A;< was given in eq. ( p.2| ).^ The expression (|3.1.1| ) applies in 
both the physical FM and AFM phases, and may be analytically continued throughout the 
respective complex-temperature extensions of these phases. 

In the physical PM phase, a similar calculation starting from (p.6| ) or ( p.lCI| ) yields 



k^'K-^C = v-Hl-v^Y/^ 



2\ 



(3.1.2) 



where v = tanhi^, and A;> was given in ( p.3[ ). Having these exact expressions for the specific 
heat, we proceed to work out its behaviour at complex-temperature singularities. 



^In passing, we note that Houtappel's expressions for the internal energy and specific heat, eqs. (108) 
and (109), respectively, in Ref. are incorrect if one uses the integrals ei(/3) and e2(/3) as he defines them, 
with the range of integration from = to 2it. If, instead, one takes the range of integration from = 
to (f) — 7r/2, so that the integrals are just the usual elliptic integrals K{yf^) and E{y/P), then his eqs. (108) 
and (109) become correct. 
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3.2 Vicinity of z = -1 

As one approaches the point z = —1 (denoted zg as above) from within either the complex- 
temperature FM or AFM phase, the specific heat diverges, with the dominant divergence 
arising from the first term in ( |3.1.1| ), which becomes —2(1 + z)^^. (There is also a weaker. 



logarithmic divergence from the term involving K{k^).) Hence, we find 

a' =a' =2 (3.2.1) 

e,FM I, AFM ^ ' 

Now, K = —(1/2) In 2;, so choosing the branch cut for the complex logarithm to lie along the 
negative real axis and choosing the first Riemann sheet for the evaluation of the logarithm, as 
z approaches —1 from above or below the negative real axis, one has Kg = ^i-K /2 respectively, 
and hence in both cases 

ks^C ^ ^7^^^ , as z^-l (3.2.2) 



It is interesting to relate the critical exponent (|3.2.1|) to the critical exponent a'^ for the 



specific heat on the triangular lattice at the point u = Ug = —1/3, which corresponds, via 
( p.l2| ), to z' = z = —1 on the honeycomb lattice. (Recall that the point m = — 1/3 is an 
endpoint of a singular line segment protruding into the complex-temperature FM phase; 
hence, it can only be approached from within this phase, and so a'^ = a'ePM-) R-^f. [[T0| , 
using the exact expression for the free energy on the triangular lattice, we calculated that 
Qfg = 1. Given the star-triangle relations which connect the Ising model on these two lattices 
and the fact that the Taylor series expansion of u + 1/3, as a function of z', in the vicinity 
of 2;' = — 1 {= z on the honeycomb lattice), starts with the quadratic term, 

« + 1 = i(l + zr + ^(1 + zf + 0((1 + zT) (3.2.3) 

it follows that the exponents = ct^^^j,^ = 2 at 2; = —1 on the honeycomb lattice 

have twice the value of = 1 at m = — 1/3 on the triangular lattice. Note, however, that 
the leading divergence in C arises from different terms on the two different lattices (in the 
triangular lattice, the leading diverg ence arises from the term proportional to E{k^^ as 
discussed in Ref. [p^ ). 

3.3 Vicinity of z = ztz 

The points z = ±i can be approached from within the complex-temperature extensions of 
the FM, AFM, and PM phases. For the approach to 2; = ±i from within the complex FM and 
AFM phases, we find from (|3.1.1| ) that the first term and the term involving E{k<) yield finite 
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contributions, while the term involving K{k^) diverges logarithmically, as ±(4i/7r)i^(fc< 
-1). Using the fact that as A ^ ±1, K{\) (1/2) ln(16/(l - A^)), and the Taylor series 
expansion of /c^ in the neighborhood oi z = ±i, 

kl = l-2{zTi? + 0{{zTiy) (3.3.1) 

we can express the most singular term on the RHS of eq. ( p.l.l|) as =F(22/7r) ln[(2; =i= i)^]. 
Evaluating K = —(1/2) In 2; for z = ±i on the first Riemann sheet of the logarithm, we have 
K = =i=i7r/4, so that 

(3.3.2) 



fco^C- ±— In 



Stt 



It follows that for z = Zs,± = ±i, 

«1,FM = (^'s,AFM = (log- div.) (3.3.3) 

These results are the same as we found for the specific heat exponent on the triangular lattice 
at the point u = —1 corresponding, via ( |2.12| ) with z' = z, to z = ±i on the honeycomb 
lattice. 

For the approach to the points v = Vs± = ±« from within the complex-temperature PM 
phase, we use the expression for C in this phase, eq. ( |3.1.2| ). We find that the term involving 
K{k^) produces a logarithmically divergence in C, so that the specific heat exponent as,pM = 
as is 

tts = {log. div.) (3.3.4) 

Taking the branch cuts for the factor (1 + 3^^)^/^ to lie along the semi- infinite line segments 
from ±i/^/3 to ±ioo, and taking the approach such that (—1)^/^ is evaluated as +i, we 
find that this term yields (2i/7r) ln[(l — k"^)]. Using the Taylor series expansion of k^, as a 
function of v, near v = i, 

kl = l-2i{v-if + 0{{v-if) (3.3.5) 
and its complex conjugate for v —>■ — z, and the result K = arctanh(±z) = ±z7r/4, we find 

777 

kn'C r^--HivT^f] (3.3.6) 

o 

(In the evaluation of the function arctanh(C) = (1/2) ln[(l + — ()] here and below, we 
again use the first Riemann sheet of the logarithm.) 
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3.4 Vicinity of v = ±z(3) 



-1/2 



We next determine the singularities of the specific heat as one approaches the endpoints 
V = ±f e = ±i/ ^/S of the semi- infinite hne segments protruding into the complex-temperature 
PM phase. Note that all directions of approach to f = ±i/^/3 except exactly (down or up, 
respectively) along the singular line segments occur from within the complex PM phase. 
From ( |3.1.2| ) we find that C is divergent, with the leading divergence arising from the term 



involving E(ky). This term gives ±{A\/3/n){l + 3^;^) ^ as f ^ ±z, so 

tte = 1 (3.4.1) 

Using K = arctanh(±z/-\/3) = ±m/6, we have 



^ ^ 33/^(1 + 3.^) ^ ^ ±7! ^'-'-'^ 

3.5 Elsewhere on the Complex- Temperature Phase Boundary 

The free energy f{K) is non-analytic across the complex-temperature phase boundaries, 
and hence, of course, this is also true of its derivatives with respect to K, in particular the 
internal energy U and the specific heat C. As an illustration, consider moving along a ray 
outward from the origin of the z plane defined by z = re*^ with 9 < 7r/2. For a given 9, as 
r exceeds the critical value rc{9), one passes from the complex-temperature FM phase into 
the complex-temperature PM phase. At the phase boundary, as is clear from Fig. 1(b), the 
elliptic modulus /c< has magnitude unity and can be written A;< = e*"^, where the angle 
depends on 9. As we discussed in connection with Fig. 1(h), z = Zc is mapped to /c< = 1, and 
z = i to = —1; (f) increases from at = to tt at = 7r/2. Hence, for < < 7r/2, 
has a nonzero imaginary part. Now when one passes through the FM-PM phase boundary 
along the ray at this angle 9, one changes the argument of the elliptic integrals from fc< = e'^'^ 
to k^ = 1/k^ = e"**^. The elliptic integrals K{k) and E{k) are analytic functions of k"^ with, 
respectively, a logarithmically divergent and a finite branch point singularity at A;^ = 1 and 
associated branch cuts which may be taken to lie along the positive real axis in the k"^ plane. 
In particular, K{k) and E{k) are both analytic at the point k = k^ = e*"^ for < 9 < tt/2. 
Hence, when we replace the argument /c< by A;>, which is the complex conjugate of A;< 
on the unit circle, we have F{ky = e"**^) = F{k^ = e^"^)* for F = K,E. Since these 
elliptic integrals are complex for generic complex k^, it follows that their imaginary part is 
discontinuous across the FM-PM boundary. The coefficients of the elliptic integrals are also 
different functions in the FM and PM phases, and these coefficients are discontinuous as one 
crosses the boundary between these phases on the above ray. Combining these, we find that 
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the specific heat itself is discontinuous as one moves across the FM-PM boundary on this 
ray. 



4 Complex- Temperature Behaviour of the Spontaneous 
Magnetisation 

The spontaneous magnetisation is given in the physical FM phase by [T^ 



;i + ^2)3/8^_4^ + ^2)l/8 ^^^^ 



(l-z)3/4(l + z)l/4 

and vanishes identically elsewhere. Observe that {1 — Az + z^) = (1 — z/ Z(^{1 — Zcz). The 
expression ( [4.1D for M can be analytically continued throughout the complex-temperature 
extension of the physical FM phase. It evidently vanishes continuously at both the physical 
critical point z = Zc, with critical exponent /3 = 1/8, and at the complex-temperature points 
z = Zs± = ±i, with exponent 

/3s = l (4.2) 

As we have observed earlier this exponent is the same as the exponent jSg^t = 3/8 
characterising the zero in the magnetisation on the triangular lattice at the point u = —1 
corresponding, via eq. (|2.12|) , to z' = z = ±i on the honeycomb lattice; however, Ps,hc = Ps,t 
differ from jSg^sq = 1/4 characterising the zero in M on the square lattice at m = — 1. This is 
a violation of universality, since this critical exponent is evidently lattice-dependent. A way 
of understanding the origin of this violation was discussed in Ref. [§] and directly reflects the 
fact that the Hamiltonian and hence the internal energy are not, in general, real numbers at 
complex-temperature singular points. 

Furthermore, M has a divergence ai z = z^ = —1, with exponent 

(3i = -\ (4.3) 

Note that the apparent zero at z = l/z^ and the apparent divergence at 2; = 1 do not 
actually occur, since these points are outside of the complex-temperature extensions of the 
FM phase, in which the (analytic continuation of the) formula ( [4.1| ) applies. We recall 
that, as a consequence of the star-triangle relation which connects the Ising model on the 



triangular and honeycomb lattices ET 



Mtiu) = Mhc{z') (4.4) 
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u, t 


z, he 


M 


(3, t 


P, he 


Uc = 1/3 


Zc = 2-y/3 


(cont.) 


1/8 


1/8 


-1 


±i 


(cont.) 


3/8 


3/8 


-1/3 


-1 


div. 


-1/8 


-1/4 



Table 1: Comparative singularities of M on the triangular (t) and honeycomb (he) lattices. 
Points are related according to the transformation ( |2.12| ) with z' = z (he). The notations 
(cont.) and div. denote, respectively, a point where M vanishes continuously and where it 
diverges. 



where g2| 



and z' (= z on the honeycomb lattice) was given in eq. ( |2.12|) . 

In Table 2 we list a comparison of singularities in Mt and Mhc- Note that, aside from the 
points listed in this table, these functions vanish discontinuously along the borders of the 
respective complex-temperature FM phases. The fact that the exponent /?£ with which Mhc 
diverges at the point 2; = — 1 on the honeycomb lattice is twice the exponent with which 
diverges at the point u = — 1/3 on the triangular lattice (which corresponds to z' = z = —1 
via (|2.12|) ) follows from the star-triangle relation connecting the Ising model on these two 



lattices and the property that the Taylor series expansion of m + 1/3, as a function of z', in 
the vicinity of 2;' = — 1, starts with the quadratic term, as given in eq. ( 3.2.3| ). 



Since the honeycomb lattice is loose-packed, one immediately infers the staggered mag- 
netisation Mhc,st from the (uniform) magnetisation Mh^'- formally, 

MhcAy) = Mhciz y) (4.6) 

Hence, Mhc,st has continuous zeros at both the physical critical point y = Uc and the complex- 
temperature point y = —1, and diverges at the points y = ±i. 



5 Analysis of the Low- Temperature Series for x 
5.1 General 

Here we shall study the complex-temperature singularities of the susceptibility x which occur 
as one approaches the boundary of the (complex-temperature extension of the) FM phase 
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from within this phase. In the next section we shall carry out a similar analysis for the 
staggered susceptibility, x^"'\ particular, we consider the behaviour in the vicinity of the 
points z = —1 and z = ±i. For the study of "we use the low-temperature series expansion 
for X, which is given by 

oo 

X = 4^3(1 + ^c„z") (5.1.1) 

ra=l 

For the study of x^^\ we use the corresponding low-temperature series expansion 



oo 



CXJ 

%3(l + ^c„,„y") (5.1.2) 



n=l 



where y = I/2; is the expansion variable in the AFM phase. These expansions have finite radii 
of convergence and, by analytic continuation from the physical low-temperature intervals 
< z < Zc and < y < yc, apply throughout the complex-temperature extensions of 
the FM and AFM phases, respectively. (Here yc is the critical point separating the PM 
and AFM phases. As can be seen from Fig. 1(a), this occurs at z = 1/zc = 2 + y^, so 
that yc = 2 — y^, the same numerical value as the critical point Zc separating the PM 
and FM phases.) Since the z^ and y^ factors are known exactly in the respective series 
expansions ( p.l.l| ) and ( p. 1.2 ), it is convenient to study the remaining factors Xr = x/(4-z^) 



and Xr"^ = X*'"V (4y'^). For the honeycomb lattice, the expansion coefficients c„ and Cn,a were 
calculated by the King's College group to order n = 13 (i.e., x ^^(^ X^""^ to 0{z^^)) in 1971 
15| , and to order n = 18 (i.e., x and x*-"^ to 0(z^^)) in 1975 |]T6|. We have checked and found 



that apparently these series have not been calculated to higher order subsequently [|T7|, |18|. 
We have analysed these series using dlog Pade and differential approximants (a recent review 
of the methods is given in Ref. |]19[)- For the approach to z = — 1 (denoted Zi) from within 
the complex-temperature extension of the FM phase, we write the leading singularity of x 

as 

Xiz) ~ A^|l - z/ze\-^'^il + ai,,|l - z/ze\ + ...) (5.1.3) 

where A'^ and 7^ denote, respectively, the critical amplitude and critical exponent, and the 
dots ... represent analytic confluent corrections. Similarly, for the approach to this point, 
which in this context we denote yi = 1/zi = —1, from within the complex-temperature 
extension of the AFM phase, we write the leading singularity of x^"'' as 

x^'Ky) ~ 4ali - y/vil'M'^ + «i/,a|i - y/vel + •••) (5.1.4) 

(As is clear from Fig. 1(a), this point z = 1/y = —1 can only be approached from within the 
complex-temperature FM or AFM phases.) As in our earlier studies of complex-temperature 
singularities of x in the Ising model on the square and triangular lattices f^, |10|, we have 
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not included non-analytic confluent corrections to the scaling form in eq. ( ^.1.3|) since, as 



we have discussed before, previous work has indicated that they are very weak or absent for 
the usual critical point of the 2D Ising model. 

Before proceeding, we consider the implications of the exact relation (|2.11|) . Given that 
Xt{u) has a singularity at m = Me = —1/3, it follows that the sum Xhd^') + xi^H-^') the 
honeycomb lattice has the same singularity at the point z' = 1 corresponding via ( p.l2|) 



to u = —1/3. But this does not, by itself, determine the singularities in the individual 
functions Xhc xi""* cit this point. If one could prove that both Xhc xi""* necessarily 
have the same singularity at = — 1, then the relation (|2.11|) , together with the result that 



Xt{u) ~ (1 + 3m) with 7g = 5/4 0, |T0|, and the Taylor series expansion of m + 1/3 as a 
function of z' (i.e., z on the honeycomb lattice) in the vicinity of 2;' = — 1, eq. (|3.2.3| ), would 



imply that 7^ = 7^,^ = 27^ = 5/2. However, although it is plausible that Xhc xi""* do 
have the same singularities at z = — 1, there is no simple relationship between the respective 



low-temperature series for these two functions, as is clear from the first few terms [fLSf , 

X = Az^ [1 + 6z + 21z^ + 122^2 + 516/ + 2148/ + ...] (5.1.5) 

and 

^(a) = 4^3j^ ^ Q . ^ ^ 3^2 ^ 2?/3 + I2y^ + 2Ay^ + ...] (5.1.6) 
Hence, an explicit series analysis is worthwhile to obtain the critical exponents. 

5.2 Exponent at z = —1 Singularity 

Since, as expected, we obtained more precise results from the differential approximants 
than the dlog Pade approximants, we shall concentrate on the former here. Our notation 
follows that of Guttmann |jl9[ and our earlier papers |10[, so we shall only describe it 



briefiy. In this method, the function / = Xr(C) being approximated satisfies a linear ordinary 
differential equation (ODE) of K'th order, £M,L/if(C) = EjLo <5i(C)^VK(C) = ^(C), where 
Qi(C) = Yle=oQj/C^ ^iid -R(C) = Y,i=oRtC (and C denotes a generic complex variable). We 
shall use the implementation of the method in which D = (d/d(. The solution to this ODE, 
with the initial condition /(O) = 1, is the resultant approximant, labelled as [L/Mq] Mk]- 
The general solution of the ODE has the form fj{C) ~ Aj{C)\C - Q\-p^ + B{C) for C Q- 
The singular points (j are determined as the zeroes of Qk{C) and are regular singular points 
of the ODE, and the exponents are given by —pj = K — 1 — QK-i{Cj)/{CjQ'KiCj))- We 
recall that if the series for the function is calculated to order , then one can compute the 
differential approximants [L/Mq] Mi] up to order L + Mq + Mi = N — 2. We shall list the 
results from the differential approximants from to Xr and xi""^ from L + Mq + Mi = 13 to the 
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[L/Mo;M^] 




\^sing 


„ / 

It 


[1/6 


6] 


-1.0096461 


9.6 X 10~^ 


2.5188 


[1/6 


7] 


-1.0009020 


9.0 X 10"^ 


2.3989 


[1/7 


5] 


-1.0056732 


5.7 X 10"^ 


2.4418 


[2/6 


6] 


-1.0047536 


4.8 X 10"^ 


2.4588 


[2/7 


5] 


-0.9945961 


5.4 X 10"^ 


2.2967 


[3/6 


7] 


-0.9959243 


4.1 X 10-=^ 


2.4136 


[3/7 


6] 


-0.9903576 


9.6 X lO"'"' 


2.3022 


[4/4 


6] 


-0.9948731 


5.1 X 10-'^ 


2.1301 


[4/6 


6] 


-0.9946423 


5.4 X 10-3 


2.3580 


[5/4 


5] 


-1.0007133 


7.1 X 10-^ 


2.1252 


[5/4 


6] 


-1.0064564 


6.5 X 10-3 


2.0704 



Table 2: Values of Zging and 7^ from differential approximants to low-temperature series for 
Xr{z). See text for definition of [L/Mo,Mi] approximant. We only display entries which 
satisfy the accuracy criterion \zsing — Z(\ < 10~^, where Z(^ — —1. 

maximum value possible using the series for these functions calculated to 0{z^^) and 0{y^^), 
viz., L + Mo + Ml = 16. We take K = 1, which will be adequate for our purposes, and use 
unbiased differential approximants since this allows us to apply an extrapolation technique 
as in our earlier work. In this technique, we plot the value of the exponent obtained from 
each differential approximant as a function of the distance of the corresponding pole location 
from the inferred exact position of the singularity (e.g., z = —1, etc.). We then extrapolate 
to zero distance of the pole from this singularity to obtain the estimate of the exponent. Of 
course, this is essentially equivalent to using biased differential approximants. We present 
our results in Table 2. 

These results yield evidence that x has a divergent singularity at z = —1, as one ap- 
proaches this point from the complex-temperature FM phase. Since the values of the ex- 
ponent from the differential approximants show considerable scatter, it is only possible to 
extract a rather crude estimate for 7^. We obtain 

7^ = 2.4 ±0.2 (5.2.1) 

This is consistent with the following inference which we shall make for the exact value of 
this exponent: 

7; = I (5.2.2) 



14 



We note that the low-temperature series does not yield as precise a determination of this 
exponent at z = — 1 on the honeycomb lattice as the corresponding low-temperature series 
did for the susceptibility exponent on the triangular lattice at the point u = —1/3 related 
to z = z' = —1 via ( p.l2| ) |10|- The critical amplitude A'^ for the susceptibility at 2; = — 1 
will be discussed below. 



6 Analysis of the Low- Temperature Series for x^^^ 

6.1 General 

The staggered susceptibility x^""^ has a well-known divergent singularity at y = yc = 2 
with low-temperature exponent 7'''*^' = 7/4. Here we analyse the complex-temperature 
singularities of this function using the low-temperature series given above in eq. ( ^.1.2| ). 

6.2 Critical Exponent of x^'^'> at z = —1 Singularity 

Our results from the differential approximants are listed in Table 3. From these we find 
strong evidence that as one approaches the point y = 1/z = —1 from within the complex 
AFM phase, x''"^ has a divergent singularity. It is interesting that the exponent values from 
these differential approximants show less scatter than those which we found for the uniform 
susceptibility. Furthermore, there is a better correlation between the value of the exponent 
and the distance of the pole location from the inferred exact value of the singularity in Table 
3, as compared with Table 2. Using our extrapolation technique, we obtain 

7;,, = 2.50 ±0.02 (6.2.1) 
This is consistent with the following exact value, which we infer: 

7i,a = I (6.2.2) 

so that, with this inference, 

7i,a = 7i (6.2.3) 



6.3 Comment on Exponent Relations at z = —1 

Using the exact results a'lFM = 2 from eq. ( ^.2.1| ) and Pe = —1/4 from eq. ([4.3|), and our 
conclusion from series analysis that j^pM = 5/2, we find that 

«tFM + 2A + 7£,FM = 4 (6.3. 1) 
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9 Q V 1 
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Z.0U40 


[3/7 
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-0.9998642 


1.4 X 10"^ 


2.4976 


[4/6 


4] 


-1.0036469 


3.6 X 10-3 


2.6153 


[4/6 


5] 


-1.0083269 


8.3 X 10-3 


2.5168 


[4/6 


6] 


-0.9983598 


1.6 X 10-3 


2.4518 


[5/4 


6] 


-1.0096131 


9.6 X 10-3 


2.4347 


[5/5 


6] 


-1.0001620 


1.6 X 10"^ 


2.1251 


[5/6 


4] 


-0.9932140 


6.8 X 10-3 


2.4105 


[5/6; 5] 


-0.9975739 


2.4 X 10-3 


2.4589 



Tabic 3: Values of ysing and 7^ „ from differential approximants to low-temperature series for 
xi^^y). Sec text for definition of [L/Mq, Mi] approximant. We only display entries which 
satisfy the accuracy criterion \ysing — yi\ < 10-^, where y^ = —1. 
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The right-hand side of eq. ( |6.3.1|) is twice the value at physical critical points. We have 
given an explanation above of why the exponents a'^pM Pe for the singularity at 2 = — 1 
on the honeycomb lattice have twice the values of the respective exponents on the triangular 
lattice, at the point u = —1/3 which corresponds, via eq. (|2.12|) to z' = z = —1 on the 
honeycomb lattice; this followed from the star-triangle relation connecting the Ising model 
on these two lattices together with the fact that the Taylor series expansion of u + 1/3, as 
a function of z', starts at quadratic order. We have also noted above the connection of our 
finding from the series analysis that 7^ for the honeycomb lattice has twice the value of the 
corresponding 7' exponent for the singularity at m = —1/3 on the triangular lattice with the 
exact relation ( p.ll|) . Since each of the exponents on the left-hand side of eq. ( |6.3.1|) has 



twice the value of the respective exponent for the corresponding singularity at m = —1/3 on 
the triangular lattice, and since this is a linear equation, the right-hand side is also twice the 
value of 2 which holds for the triangular lattice. 

One may also consider the analogous equation for the approach to z = —1 from within 
the complex-temperature extension of the AFM phase. We have extracted the exact value 
a'i AFM = 2 in eq. ( |3.2.1|) and, as discussed above, given the loose-packed nature of the 



honeycomb lattice and the resultant relation ([4.6|) , it follows that the staggered magnetisation 
diverges with the exponent Pi^st = /?£ = — 1/4 as one approaches the point z = —1 from within 
the complex-temperature AFM phase. Combining these exact results with the conclusion 
from our analysis of the low-temperature series for the staggered susceptibility that ^ = 
5/2, we find 

+ 2A,,t + 7L =4 (6.3.2) 



a, 



,AFM 



in complete analogy with eq. (|6.3.1|) , as expected for a loose-packed lattice. 



6.4 Critical Amplitude for x^"^ at z = —1 

In order to calculate the critical amplitude A'^^ in the staggered susceptibility as one ap- 
proaches z = y = —1 from the complex AFM phase, we compute the series for {x^f^Y^'^'^''' 
using our inferred value 7^^ = 5/2. Since the exact function [x^r^Y^'^^''' has a simple 
pole at z = —1, one performs the Fade analysis on the series itself instead of its loga- 
rithmic derivative. The residue at this pole is —yi{A'^^^Y^"''^'°-, where A'^ar denotes the 
critical amphtude for x^"^ Extracting A'^^r multiplying by the prefactor, we obtain 

a; ^ = -0.700 ±0.010 (6.4.1) 

An analysis of the series for xl-^^'^ to extract the critical amplitude for the (uniform) sus- 
ceptibihty as one approaches z = —1 from within the complex FM phase did not yield precise 
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results, presumably because of the shortness of the series. However, having already inferred 
that X and x''"'' have the same power- law divergence at z = ?/ = — 1, as approached from the 
complex FM and AFM phases, respectively, we can use the relation (|2.11|) to compute A'g^ 
indirectly. For this purpose, we recall that on the triangular lattice, at the corresponding 
point M = Me = — 1/3, the (uniform) susceptibility Xt has the leading singularity P, [l^ 



A' (l + 3n) 



-5/4 



(6.4.2) 



Using this, together with the Taylor series expansion ( p. 2. 3 ), we find the following relations 
among the critical amplitude j at m = — 1/3 on the triangular lattice and and A'g^^ at 
2; = — 1 on the honeycomb lattice: 



2 3^/^ A 



e,i 



We next use our previous determination of A'^ ^ 



A' 



e,t 



[0 

0.05766 ±0.00015 



(6.4.3) 



(6.4.4) 



which agrees with, and has somewhat smaller uncertainty than, an earlier determination of 
this quantity by Guttmann 0. Substituting (|6.4.1|) and ( |6.4.4|) into ( p.ll| ), we obtain the 
critical amplitude for the uniform susceptibility. 



A[ = 0.245 ± 0.010 



(6.4.5) 



7 Other Singularities 

We also used the low-temperature series for x and x^""^ to investigate the singular behaviour 
of these functions as one approaches the points z = Zs± = ±i from within the complex- 
temperature FM and AFM phases, respectively. However, we were not able to obtain con- 
clusive results. This is similar to our previous experience investigating the behaviour of x on 
the triangular lattice in the vicinity of the singular point m = — 1 as approached from within 
the complex FM phase. These points, i.e., z = ±i on the honeycomb lattice and -u = — 1 or 
equivalently z = ±i on the triangular lattice, share in common the property that they are 
intersection points where two arcs of the complex-temperature phase boundary curves cross. 

Concerning singularities of x and x^""^ elsewhere on the boundaries of the complex- 
temperature FM and AFM phases, respectively, it is quite possible that these functions 
may exhibit discontinuities, as we found for C and M. As is well-known from work on first- 
order phase transitions, the analysis of the low-temperature series by itself, and similarly, 
the high-temperature series by itself, does not, in general, provide a sensitive probe for such 
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discontinuities. Of course if one had sufficiently long series, the comparison of the limits as 
one approached a phase boundary from, say, the FM and the PM phases using the respective 
low- and high-temperature series, could be of value in this regard. 



8 Behaviour of % in the Symmetric Phase 

The theorem proved in Ref. and discussed further in Ref. implies that, for the Ising 
model on the square lattice, x has at most finite non-analyticities as one approaches the 
boundary of the complex-temperature extension of the PM phase, aside from the physical 
critical point aX v = Vc- We would expect a similar theorem to hold for the honeycomb 
lattice, although to show this with complete rigour, it would be desirable to perform an 
analysis of the asymptotic behaviour of the general connected 2-spin correlation function 
< o"o,oO"m,n > as r = (m? + 71^)^^"^ oo for this lattice, which has not, to our knowledge, 
been done. Assuming that such a theorem does hold, it would follow, in particular, that, 
x{v) and x^"'\v) = x{~v) would have finite non-analyticities as one approaches the points 
V = ±i from within the PM phase. We have analysed the high-temperature series expansion 
for to investigate the singularities at f = ±i. This series is of the form 

oo 

X=l + E«n^" (8-1) 

n=l 

It has a finite radius of convergence and, by analytic continuation from the physical high- 
temperature interval < v < Vc, applies throughout the complex extension of the PM 
phase. The high-temperature series expansion ( ^.11) is known to 0(f^^) [^, Since we 



anticipated a finite singularity, we analysed this series using differential approximants, which 
are capable of representing this type of singularity in the presence of an analytic background 
term. However, as we recall from Fig. 1(c), there are two semi-infinite line segments which 
protrude into the complex-temperature PM phase, with endpoints at f = ±f g = ±i/ Vs. We 
found that the series are not sensitive to the singularities at v = ±i, presumably because of 
the effect of the intervening singular line segments and their endpoints at f = ±i^/3. We 
have also tried to study the singularities in x at these endpoints. Again, our study did not 
yield an accurate value for the exponent 7e, presumably due to the insufficient length of 
the series. However, the (scattered) values of 7e were consistent with the expectation that 
7e < 0. 
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9 Conclusions 



In this paper wc have investigated complex-temperature singularities in the Ising model on 
the honeycomb lattice. As part of this, we have discussed the complex-temperature phases 
and their boundaries. From exact results, we have determined these singularities completely 
for the specific heat and the uniform and staggered magnetisation. Prom an analysis of low- 
temperature series expansions, we have found evidence that x and x^"^ both have divergent 
singularities z — —1 = Z(, (where z — e"^^), with exponents 7^ = „ = 5/2. The critical 
amplitudes at this singularity were calculated. We have found that the linear combination 
of exponents a' -\-2l3 -\- 7' is equal to 4 rather than 2 at 2; = —1. The connection of these 
results to corresponding complex-temperature singularities on the triangular lattice was dis- 
cussed. Finally, some results on complex-temperature singularities reached from within the 
symmetric phase were given. 

One of us (RS) would like to thank Profs. David Gaunt and Tony Guttmann for infor- 
mation about the current status of the series expansions for the honeycomb lattice. This 
research was supported in part by the NSF grant PHY-93-09888. 
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Figure Caption 



Fig. 1. Complex-temperature phases and associated boundaries in the variables (a) 
z, (b) A;<, and (c) v. FM, AFM, and PM denote ferromagnetic, antiferromagnetic, and 
paramagnetic, Z2-symmetric phases. In Fig. 1(b), the semi- infinite line segment extends 
along the negative real axis from /c< = to /c< = — oo. In Fig. 1(c), the semi-infinite line 
segments extend along the imaginary axis from v — ±i/-\/3 to ±ioo, respectively. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2076v 1 
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